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ABSTRACT 

Excess soft X-ray emission in clusters of galaxies has so far been detected for 
sources that lie along lines-of-sight to very low Galactic HI column density (such 
as Coma, A1795, A2199 and Virgo, Nh ~ 0.9-2.0 x 10^° cm'^). We show that the 
cluster soft excess emission can be investigated even at higher Nh, which provides 
an opportunity for investigating soft X~ray emission characteristics among a large 
number of clusters. 

The ROSAT PSPC analysis of some members of the Shapley concentration 
(A3571, A3558, A3560 and A3562, at iV^/ ~ 4 - 4.5 x lO^o cm'^) bears evidence 
for excess emission in the 1/4 keV band. We were able to confirm the finding 
for the case of A3571 by a pointed SAX observation. Within the current sample 
the soft X-ray flux is again found to be consistently above the level expected 
from a hot virialized plasma. The data quality is however insufficient to enable 
a discrimination between alternative models of the excess low energy flux. 

Subject headings: galaxies:clusters:individual (A3571, A3558, A3560, A3562) - 
intergalactic medium 

1. Introduction 

Clusters of galaxies are strong emitters of soft X-rays, often in excess of the contribution 
from the hot intra-cluster medium (ICM; Lieu at al. 1196a; Lieu et al. 1996b; Mittaz, 
Lieu and Lockman 1998; Lieu, Bonamente and Mittaz 1999; Bonamente, Lieu and Mittaz 
2001). The PSPC detector onboard ROSAT is well suited for the observation of diffuse 
sources in the 1/4 keV band: along with a large effective area (of about 200 cm^ at 0.28 
keV) , the large field-of-view of 1 degree radius allowed the contemporaneous acquisition of 
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source signals (usually near boresight) and local background (at large detector radii) with 
equal exposures. The Low-Energy Concentrator Spectrometer (LEGS) aboard BcppoSAX 
also has appreciable effective areas (of ~ 10 cm^ at 0.28 keV) and can be used for similar 
investigations. 

Soft X-ray fluxes are significantly absorbed by intervening Galactic material (mainly 
HI and Hel) along the line of sight. At the Galactic Nh values of the Shapley supercluster 
(rv^ 4 — 4.5 X 10^° cm~^), 0.2 keV photons encounter a total optical depth of ~ 4, which 
accordingly reduces their flux by about 98 %. The Hel absorption cross-section was recently 
revised (Yan et al. 1998 and references therein), leading to a reinstatement of the original 
absorption code of Morrison and McGammon (1983) as the most rehable. This code is 
therefore employed in the present work. The crucial task, however, is to measure Nh with 
a spatial resolution < the cluster size. For details on how this was fulfilled, see section 2. 
In this paper we present results on a sample of nearby clusters (A3571, A3558, A3560 and 
A3562, at redshift z~ 0.04-0.05) in the Shapley concentration. The analysis of four PSPG 
observations revealed statistically significant 1/4 keV band (~ 0.2-0.4 keV) excess emission 
for all of the clusters, although the intrinsic low surface brightness of A3560 and A3562 
imphed a hmited S/N. Data from two pointed SAX/LEGS observations were available for 
two of the four clusters (A3571 and A3562); clear evidence of soft excess emission is present 
in the A3571 data, while results for A3562 were marginal, due to the faintness of the 
cluster. Thus the soft excess syndrome may be common to a large fraction of clusters, as 
suggested also by our earher analyses of a nearby sample of clusters with EUVE, contrary 
to the findings of Bowycr ct al. (1999,2001) which claim no general importance of the effect 
(see also Lieu et al. 1999 and Bonamente, Lieu and Mittaz 2001 for reinforcement of the 
original conclusions through confirmatory re-observations conducted in such modes as to 
address all the objections of our critics.). This paper presents the data analysis method and 
results for the current sample, with emphasis on the PSPG analysis; future works will be 
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devoted to larger scale surveys, with both PSPC and LEGS instruments, in order to assess 
the cosmological impact of the phenomenon. 

2. Infrared, 21— cm and X— ray observations of A3571, A3558, A3562 and 

A3560 

The Shapley Supercluster is a large region of galaxy overdensity around a — 13'*25"*, 
5 = -30° and redshift ~ 0.045 (Shapley 1930; Ettori, Fabian and White 1997 and 
references therein) comprising more than 20 galaxy clusters. We focused our attention 
on four Supercluster members: A3571 and A3558 are rich and hot clusters (Markevitch 
et al. 1998; Markevitch and Vikhlinin 1997), A3560 (Ebeling et al. 1996) and A3562 
(White 2000; see also Ettori et al. 2000) are within ~ 2 degrees from the former and show 
lower X-ray luminosity (Fig. 1). All clusters were targets of pointed PSPC observations 
(archival identification code RP800287 for A3571, exp. time 5 ks; RP800076 for A3558, 28 
ks; RP800381 for A3560, 5 ks; RP800237 for A3562, 18 ks), and two of them of pointed 
SAX LECS/MECS observation (observation number 60843002 for A3571, 29 ks of LEGS 
exposure, and 60638001 for A3562, 19 ks) . 

Dedicated narrow-beam 21-cm observations of galaxy clusters are available from the 
NRAO 100-meter GBT telescope ^ with a resolution of 10 arcmin (Murphy, Sembach 
and Lockman 2001; see Table 1). In addition we also consulted the Dickey and Lockman 
(1990) survey data at lower resolution and IRAS 100 //m maps with ~ 3 arcmin resolution 
(Wheelock et al. 1994). In particular, the IRAS data has a two-fold relevance to soft 
X-ray studies: on one hand the positive N^^-IR correlation allows the use of 100 /xm maps 



^The National Radio Observatory is operated by Associate Universities, Inc., under co- 
operative agreement with the National Science Foundation. 
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as tracer of the total column density of Galactic material (Boulanger and Perault 1988); 
on the other hand, the Nj^-soft X-ray background (SXB) anticorrelation (Snowden et al. 
1998) requires precise knowledge of the HI column density at the position where the local 
background was determined. In the case of PSPC this is typically a few x 10 arcmin away 
from the cluster center. 



3. PSPC analysis 

The analysis of PSPC data was performed following the prescriptions of Snowden et al. 
(1994); in particular, detector gain corrections were applied using the dedicated FTOOL 
PCPICOR, and detector vignetting was accounted for in the background subtraction 
procedure. The spectral analysis involved only PI channels 20-201 (~ 0.2-2.0 keV by 
photon energy), in order to avoid any calibration uncertainties of the lower energy 1-19 PI 
channels. Following the nomenclature reported in Snowden et al. (1994), PI channels 20-41 
will be referred to as the R2 band (also 1/4 keV or C band), and channels 42-69 as R34 
band. 

First, we compared the IR emission of the cluster regions with the emission in the 
corresponding background regions (Fig. 1); deviations were always within at most 10 %, 
and this implies that the SXB of the background regions is indeed relevant to the cluster 
regions without further scaling (see Snowden et al. 1998). Secondly, a comparison of the 
IR emission between the cluster regions covered by the GBT beam (~ 10 arcmin angular 
diameter) and the outer cluster regions showed that the Nh is smoothly distributed on the 
cluster's scale. In accordance to the N^^-IR correlation (with a scaling factor of about 1.2 
xlO^° cm~^ (MJy sr~^)~^; Boulanger and Perault 1988) the 100 iim maps imphed only 
marginal adjustments from the central measurements; as a consequence, we decided to 
apply Galactic column densities as in Table 1 for all subsequent analyses. 
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The limited passband of the PSPC instrument can only be used to measure approximate 
values of the cluster hot gas temperatures. Thus we modelled the A3571 and A3558 
spectra for annular regions with the best-fit temperatures of ASCA (Markevitch et al. 
1998; Markevitch and Vikhlinin 1997; in agreement with SAX/MECS results for A3571, 
see Nevalainen et al. 2001), with the exception of the central regions where temperature 
gradients are too steep for the point spread function (PSF) of ASCA: here T was fitted 
to the data. Temperatures of A3560 and A3562 were determined by spectral fits to the 
PSPC data, and results agreed with previously pubhshed measurements (White 2000 and 
Ebeling et al. 1996). We modelled the entire PSPC X-ray spectrum (0.2-2.0 keV) with a 
one-temperature photoabsorbed thin-plasma emission code (WABS and MEKAL codes in 
XSPEC). 

Background spectra for subtraction purposes were extracted from same pointed 
observation as the cluster's, and from annular regions immediately outside the detected 
cluster emission (see Fig. 1); statistical errors in the background were properly propagated. 
This provides a spatially and temporally contiguous (or in situ) background. Best parameter 
values for the hot gas are shown in Table 1. As the values indicate, some of the fits 
are poor, due to the excess emission above the 1-T model at low energies. In more detail, 
we compared the best-fit model with the measured emission in the R2 and R34 bands 
(Fig. 2 and 3). The 1/4 keV band results show consistently positive residuals for all of 
the 4 clusters; the R34 band bears, on the other hand, evidence of depleted fluxes in some 
regions. For A3571 and A3558, there is evidence of an increasing relative importance of 
the excess component with cluster radius, an effect already found in EUVE data of other 
galaxy clusters (Virgo, A1795 and A2199; see, e.g., Bonamente, Lieu and Mittaz 2001 and 
Lieu, Bonamente and Mittaz 1999). 
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4. SAX LEGS analysis 

The SAX mission lends an additional view of the phenomenon, with an instrument, 
LEGS, that employs different and complementary characteristics. With a field of view of 
radius ~ 16 arcmin, local background can be obtained only for compact clusters (such 
as A3562); a long-term campaign of accumulation of high Galactic latitude 'blank' fields, 
however, enables one to establish a template detector background applicable to clusters 
which fill the entire FOV (such as A3571; see Parmar et al. 1999). The energy dependent 
Point Spread Function (PSF) of the LEGS instruments complicates the analysis of extended 
objects such as clusters of galaxies, as its 90% encircled energy radius increases from 6' 
to 9' between 8 keV and 0.2 keV. In this work, we accounted for the PSF by generating 
appropriate instrument response files to be used when fitting the accumulated spectra. For 
the response generation we used the ray-tracing code LEMAT (Lammers 1997), taking into 
account the extraction region size and source position on the detector, mirror vignetting 
and obscuration by the detector support structure (strongback) . ^ 

Modeling of the low-energy emission from A3571 was performed through a joint fit of 

MECS/LECS spectra to a photoabsorbed MEKAL code (as for PSPC), which returned the 
best-fit parameters of the hot IGM consistent with those of Table 2 (see Nevalainen et al. 

^For any given spatial region and energy range, LEMAT estimates the contribution of 
photons from regions external to the extraction region, and the number of photons originating 
from within the extraction region that are detected externally (see Kaastra et al. (1999); 
Nevalainen et al. (2001) for the application of LEMAT in the case of the clusters A2199 
and A3571). Further applications of LEMAT to spectra of extended sources is in progress 
(Gosterbroek 2001); prehminary results indicate that using the LEMAT-produced responses 
in spectral fits of spectra obtained within a radius of 15' removes the significant instrumental 
distortions (except for the whole band normalization which can vary by ± 30%). 
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2001 for the details on the LECS/MECS analysis of this data set): excess emission in the 
0.15-0.3 keV region^ is again detected with high statistical significance (Fig. 4). The faint 
emission from the A3562 cluster was likewise assessed by modelling LEGS spectra with the 
customary photoabsorbed MEKAL code, where the hot ICM parameters were fixed at the 
values of Table 2; results are marginahzed by the lower S/N (see Fig. 4). It is worthwhile to 
point out that the in situ background, available for this cluster from an outer annulus where 
the source emission ceased, is fully consistent with the template 'blank' field background. 
This is a crucial test which verifies the rehabihty of the latter for those cases where the size 
of the cluster renders the former unavailable. 



5. Interpretation of R2 and R34 band fluxes 

The results shown in Fig. 2 (and Fig. 4) indicate that the 1-T model employed to 
describe the ICM falls short of the detected soft emission, revealing the soft excess emission. 
The data quality for A3571 and A3558 warrants the employment of a more complex model 
to try to interpret the detected spectral features. 

The A3558 soft excess is virtually unaccompanied by a R34 trough, and therefore a 

^At 0.28 keV, the LEGS energy resolution is ~ 32% (Parmar et al. 1997), about 3 times 
better than PSPG's. This renders it difficult to translate PSPG pulse-invariant (PI) channels 
into corresponding LEGS band for the purpose of comparison. Our present choice of using 
0.15-0.3 keV band of LEGS as equivalent to PSPG PI channels 20-41 is based on the 10% 
peak response studies of Snowden et al. (1994). Similarly, PSPG R34 band has 10 % of peak 
response 0.2-1.01 keV; moreover LEGS data of A3571 also bear evidence of depleted fluxes 
around 0.5 keV (Nevalainen et al. 2001), though no attempt is made here to extract and 
compare such fluxes with those of the PSPG. 



-9- 



second optically-thin plasma component at lower temperature can model the excess. We 
shall present the 3-6 and 6-9 arcmin annuli as representative regions, just outside the 
central cooling flow, and where signals are of sufficient statistical quality to apply detailed 
modelling, see Table 2. We likewise tested the possiblity of the excess originating as an 
inverse Compton (IC) effect, with the result that the data equally well accomodates the 
non-thermal interpretation (see Table 2 for details). We also calculated the energetic 
requirements of the additional model, with the results that the luminosity of either 
thermal/non thermal component accounts for only a fraction of the hot ICM luminosity in 
the 0.1-0.4 keV band. 

Moving to A3571, we again describe a representative region, the 4-8 arcmin annulus. 
Here the ~ 3a R2 band excess pairs with a depleted emission in the neighboring R34 band 
(Fig. 2 and 3.), i.e. the employement of a second thin plasma model at lower T brings only 
a modest improvement to the fit. We know, however, that warm photoionized gas (T ~ 10^ 
K) may present substantial opacity to soft X-ray photons (e.g., Krohk and Kallman 1984), 
and the thin-plasma approximation may not necessarily apply. A phenomenological way of 
parametrizing such behaviour is the addition of a second thin-plasma model modified by an 
absorption edge, which in fact proves a successful fit to the data (see best-fit parameters in 
Table 2). Other regions of Table 2 show a qualitatively similar behavior, see Fig. 2 and 3; 
statistical quality of the data does not however warrant a more detailed modelling for them. 

We also entertained the notion of a soft X-ray excess as pecuhar variations of the 
Galactic Nh (Arabadjis and Bregman 1999), by applying the same 1-T thermal models 
above leaving the Nh to be fitted by the data. The A3571 and A3558 PSPC spectra would 
in such case require a Galactic column density at odds with the measured values, and some 
1/4 keV band excess would nonetheless remain unexplained (see Table 1 and Fig. 5). 
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6. Discussion and conclusions 

The analysis we present in this paper shows that, at soft X-ray energies, PSPC spectra 
of the clusters in our sample (A3571, A2558, A3560 and A3562) are not well described 
by a simple one-temperature thin plasma model. Similar conclusions hold for the LEGS 
data of A3571 and A3562. The case of A3558, the cluster with the higest S/N data of 
the sample, indicates that the soft excess can be interpretated as having a thermal origin 
(from a second phase of the ICM at lower temperature) as well as the non-thermal. Earlier 
EUVE observations led to the discovery of soft excess emission for the brightest clusters 
with low Galactic HI columns (such as Goma and Virgo, see references in section 1), but 
they lack the necessary sensitivity to enable a reliable assesment of fainter sources, such as 
those reported in this work. The PSPG and LEGS data presented here show that, given a 
sufficient sensitivity, low energy spectra of galaxy clusters consistently show excess emission 
above the hot IGM contribution, in sharp contrast with the conclusions of Bowyer et al. 
(2001). Settlement of the issue regarding the origin of the emission has to await forthcoming 
observations by detectors of higher resolution. 

The BeppoSAX satellite is a joint Italian-Dutch programme. We thank the staff of the 
BeppoSAX Science Data and Operations Gontrol Genters for help with these observations. 
MB and RL acknowledge the support of NASA/ADP. JN acknowledges an ESA Research 
Fellowship. 
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Figure captions 



Figure 1: Schematic of the field of the four clusters (A3571, A3558, A3560 and A3562). 
Large sohd circles represent the approximate cluster positions, the two small circles between 
A3562 and A3558 are the groups of galaxies SC1329 and A1327 (respectively from left to 
right; see Ettori et al. 1997 for more detailed maps). Each dotted annulus (or sector thereof) 
represents the approximate region from which the background for the corresponding cluster 
was extracted; care was taken in avoiding regions occupied by SC1329 and SC1327 and 
obvious point sources. 

Figure 2: Radial distribution (in arcmin, x-axis) of the fractional R2 band excess, f], 
of the four clusters, defined as 77 = ip2 — (I2) / (I2, where p2 is observed R2 band fiux and q2 
the expected R2 band fiux from the hot ICM. Negative values of ( indicate that detected 
emission is below the expected value. 

Figure 3: As in Figure 2, except now for the fractional R34 band excess ( — (P34— Q'34) / 534 
where pu is observed R34 band fiux and ^34 is expected R34 band fiux from hot ICM. 

Figure 4: Radial distribution of the SAX/LECS fractional excess in the 0.15-0.3 keV 
band for A3571 and A3562, plotted in the same manner as Figure 2. 

Figure 5: As in Figure 2, except here the expected ICM fiuxes are derived by fitting 
Nfj. Soft excess fiuxes in some regions of A3571 remain nonetheless unexplained. In order 
to eliminate such excesses, we estimate that a Nh of 4.05 ±o!2^ x 10^° (4-8 arcmin region) 
and 3 ±q;4 x 10^° (12-16 arcmin region) would be required, significantly lower than the 
measured values (see Table 1). 
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Tabic 1: Best parameters of the single temperature hot gas for the clusters in the sample. 
Errors are 90 % confidence level (x^ + 2.7 criterion); where no errors are reported, the 
parameters were not fitted to the data (see text). Fitted T values are obtained keeping NH 
fixed at the Murphy, Sembach and Lockman (2001) values. In bold-face are the Nh values 
from pointed observations of Murphy, Sembach and Lockman (2001); the survey values are 
from Dickey and Lockman (1990), with average over a 3x3 square degrees area. 



Cluster region (arcmin) Nij Survey Nij best-fit N^j T (keV) A (d.o.f) 
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A3562 0-3 4.0 3.85 4 ±g:| 2.7 ±g:| 0.6 203.6(165) 

3-6 " " 3.56 ±1% 4.5 ±\\ 0.3 134(161) 

6-12 " " 3.65 ±%f 3.6 ±J| 0.3 205.4(163) 



A3560 0-3 4.7 4.6 5.3 ±1% 2.2 ±J| 0.38±[];i 52.7(47) 

3-6 " " 4.2 ±1% 4.2 ±f^^ 0.3 85(58) 

6-12 " " 4.1 ±11 l.Q±li 0.25±g:^l 104.4(78) 

a - Use of best-fit values kT=7.5, A=0.5 (Markevitch et al. 1998) will result in no changes to 
Fig. 2 and Fig. 3. 
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Table 2: Modelling the PSPC spectra of some regions of A3571 and A3558 with an additional 
thermal or non-thermal component; parameters of the hot ICM were fixed at the best values 
as given in Table 2 except the normalization constant (emission measure) which was fitted to 
the data. For A3558, a second MEKAL component is added to the model, with the same ele- 
mental abundances as those of the hot phase; alternatively, a power-law model is added with 
its differential photon number index fixed at the value of Q!=1.75, corresponding to emitting 
electrons having the Galactic index of cosmic rays (see, e.g.. Lieu et al. 1999). For A3571, 
the model employed was, in XSPEC language, a WABS*EDGE*(MEKAL+MEKAL). This 
model is only intended as a preliminary means of accounting for the observed spectral fea- 
tures, as a model that includes details of radiative transport and opacity is not currently 
available. These spectral models were also satisfactorily applied to the LEGS data, showing 
that PSPC and LEGS give consistent results. 
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